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ABSTRACT 

We present a simple and intuitive approximation for solving perturbation theory (PT) 
of small cosmic fluctuations. We consider only the spherically symmetric or monopole 
contribution to the PT integrals, which yields the exact result for tree-graphs (e.g. at 
leading order) . We find that the non- linear evolution in Lagrangian space is then given 
by a simple local transformation over the initial conditions, although it is not local in 
Euler space. This transformation is found to be described the spherical collapse (SC) 
dynamics, as it is the exact solution in the shearless (and therefore local) approxima- 
tion in Lagrangian space. Taking advantage of this property, it is straightforward to 
derive the one-point cumulants, ^j, for both the unsmoothed and smoothed density 
fields to arbitrary order in the perturbative regime. To leading order this reproduces, 
and provides with a simple explanation for, the exact results obtained by Bernardeau 
(1992, 1994). We then show that the SC model leads to accurate estimates for the 
next corrective terms when compared to the results derived in the exact perturbation 
theory making use of the loop calculations (Scoccimarro & Frieman 1996). The agree- 
ment is within a few per cent for the hierarchical ratios Sj = S,j/S,2~^- We compare 
our analytic results to N-body simulations, which turn out to be in very good agree- 
ment up to scales where ct ~ 1. A similar treatment is presented to estimate higher 
order corrections in the Zel'dovich approximation. These results represent a powerful 
and readily usable tool to produce analytical predictions to describe the gravitational 
clustering of large scale structure in the weakly non-linear regime. 

Key words: cosmology: large-scale structure of Universe-cosmology: theory-galaxies: 
clustering- methods:analytical- methods: numerical. 



1 INTRODUCTION 

One of the important problems in Cosmology is the origin of 
galaxies, clusters, and structures seen on huge scales in the 
spatial distribution of galaxies and in the microwave back- 
ground radiation. Given some initial seed fluctuations we 
would like to be able to make analytic predictions to de- 
scribe the gravitational clustering that generates the large- 
scale structure that we see today. 

The probability distribution function contains all the 
statistical information concerning the cosmic fields: density 
5 and velocity v fluctuations. Here we will concentrate on its 
moments, the variance, skewness, kurtosis and so on. For a 
given set of moments in the initial conditions (IC), we would 
like to derive the flnal moments using the exact dynamics 
that rule the evolution of the underlying field. The problem 
is that the exact solution to the dynamical equations is only 
known for the linear regime of the evolution of these fields 
and only approximate solutions are known for their corre- 



sponding non-linear stages. The spherical infall model has 
been used to provide a qualitative picture of how isolated 
(gravitationally bound) structures evolve up to their final 
collapse. In particular, the spherical model gives an estima- 
tion of when the first astrophysical objects formed purely 
under infiuence of gravity (see Pebbles 1980, and references 
therein). This motivated the development of the so-called 
Press-Schechter formalism (Press & Schechter 1974) which 
gives the distribution of collapsed objects as a function of 
their mass. It is commonly assumed that highly non-linear 
objects such as galaxies form at peaks {i.e., values above a 
threshold) of the underlying matter density field, what pro- 
vides a possible mechanism for the segregation (bias) be- 
tween luminous and dark matter in the universe (Kaiser 
1984, Bardeen et al. 1986). More recently, Mo & White 
(1996) have extended the spherical collapse model to build a 
model for the spatial clustering of non-linear objects such as 
dark matter halos. Here we shall explain how the spherical 
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collapse model also appears in the context of Cosmological 
Perturbation Theory. 

Cosmological Perturbation theory (PT) provides a 
framework to study small departures from linear theory: the 
weakly non-linear regime. The leading order contribution to 
the skewness for Gaussian initial conditions (GIC) was ob- 
tained by Peebles (1980). Fry (1984) extended this result to 
higher order moments and Bernardeau (1992, hereafter B92) 
found the generating function to the leading order contribu- 
tion: the tree-level solution. 

Comparison with observations and simulations made 
it necessary to develop PT for the smoothed fields. The 
smoothing corrections to the unsmoothed amplitudes for a 
power-law were first computed for the skewness for either 
a top-hat or a Gaussian window function by Juskiewicz et 
al. (1993). For a top-hat filter, Bernardeau 1994a (hereafter 
B94a) developed the machinery to systematically derive the 
smoothed hierarchical amplitudes for either the density {S.j) 
or the velocity divergence (T./) fields to an arbitrary order 
(and arbitrary power spectrum). These results are in ex- 
cellent agreement with numerical simulations {e.g., Baugh, 
Gaztaiiaga, Efstathiou 1995, hereafter BGE95; Gaztanaga 
& Baugh 1995, Colombi et al. 1996). Work by Matsubara 
(1994) and Lokas et al. (1995) include some results for the 
Gaussian-smoothed density and velocity cumulants. 

The success of the PT approach made it plausible to go 
further and study higher order (loop) corrections and the 
case of non-Gaussian initial conditions (NGIC). Despite the 
impressive achievements in the diagrammatic approach to 
loop corrections in PT by Scoccimarro & Frieman 1996a and 
1996b, (hereafter SF96a and SF96b, respectively), the ana- 
lytic entanglement faced when carrying out the calculations 
for the smoothed one-point cumulants is enormous. More- 
over, regularization techniques must be used to evaluate the 
loop corrections through the kernels in Fourier space which, 
for a value of the spectral index n — —1 lead to logarithmic 
divergences which become powers of the cut-off used in the 
regularization for n > — 1. This puts constraints on the do- 
main of applicability of the exact PT itself. When it comes 
to calculating the statistics of the smoothed fields, further 
analytic entanglement makes comparison with simulations 
and observations be restricted to few discrete values of the 
spectral index (Scoccimarro 1997, hereafter S97). For non- 
Gaussian initial conditions we have to face similar problems 
to start with, as loop corrections could enter at the same 
order than the tree-level contribution. 

The general results derived by Bernardeau (B92, B94a 
and Bernardeau 1994b, hereafter B94b) show that the gen- 
erating function follows the equation for the spherical col- 
lapse model. In this context this remarkable result lacked a 
satisfactory interpretation. What is the connection between 
PT and the spherical collapse model? Why does the gen- 
erating function follow an equation for a fluctuating field? 
Why is the smoothed generating function given by the same 
unsmoothed result but with a different argument? Is this 
related to the Gaussian nature of the IC or the flatness 
of space? In this paper we will present a simple interpre- 
tation for these results, which will allow us to extend it to 
higher-order corrections and non-Gaussian initial conditions 
(NGIC). 

This paper is organized as follows. We present the 
monopole approximation to PT in 93. In section §pl we give 



a brief account of the estimation of the cumulants in PT 
and some of its nonlinear approximations. We investigate 
the Spherical Collapse (SC) model and its relation to PT in 
iM. In section §H we present a comparison of the SC model 
with exact PT and N-body Simulations. General formulas 
for the one-point cumulants from a local transformation of 
the density field in Lagrangian space (for GIC) are provided 
in Appendix H. A comparison with analog results in Euler 
space is given in Appendix H An alternative derivation of 
the smoothing effects in Fourier space for a top-hat window 
is illustrated in Appendix |c[ Further results are derived for 
the Spherically Symmetric Zel'dovich Approximation in Ap- 
pendix O. 



2 PERTURBATION THEORY AND THE 
MONOPOLE APPROXIMATION 

We will now review briefly the main features of perturba- 
tion theory (PT) to present the spherical symmetric solu- 
tion or monopole approximation. Throughout this paper we 
shall focus on Gaussian initial conditions and the Einstein- 
deSitter scenario (FRW with f2 = 1, A = 0) to gain sim- 
plicity. However, a generalisation of the current approach to 
non-Gaussian IC and non-flat FRW universes will be pro- 
vided in two accompanying papers (see Gaztanaga & Fos- 
alba 1998, Fosalba & Gaztafiaga 1998, Papers II & III, here- 
after) . 



2.1 Equations of Motion 

In a non-relativistic presureless (coUisionless) expanding 
universe, the fleld equations for the local fluctuation field 
5{x,t) in Euler space are: 



dV a 1, „^ 

-^ + -V+-{V-\7) V 
at a a 



a 



V^ = ^TiGpa^S 



(1) 



(2) 



which correspond to the mass conservation, Euler and Pois- 
son equations respectively. The expansion factor a = a{t) = 
(1 -I- z)~^ obeys: 



H ^ Ho [no/a-' + (1 - Oo - nA)/a^ + Qa] 
8-kGpo „ A 



1/2 



3Hi ' 



3//2- 



(3) 



Here Ho and po are the present values of the Hubble constant 
and background density, and A is the cosmological constant 
(see e.g., Peebles 1980). 



2.2 Perturbation Theory Solutions 

A convenient approach to solve the field equations above is 
to expand the density contrast S assuming it is small as is 
usually done in the context of perturbation theory (PT) in 
Euler space, 

5{X, t) = 5i {X, t) + hiX, t) + 53{X,t)..., (4) 
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were we assume that Sn « 5„-i. The first term 5\ = 5i 
is the solution to the hnearized field equations and is given 
in terms of the linear growth factor D{t). The second term 
§2 ~ (5f is the second-order solution, obtained by using the 
linear solution in the source terms, and so on. 

The dominant mode to the linear PT solution, 5;, is 
given in terms of the linear growth factor D{t): 



5,(x,t) = D{t) S{x,0) 



D{t) ^ - da d~^, 

"- Jo 



(5) 



which for fl = 1, the so-called Einstein-deSitter universe, 
gives D{t) — a{t). The linear approximation is therefore triv- 
ially local: the evolution of a density fluctuation at a given 
point is not affected by its neighboring density fluctuations, 
and all fluctuations grow with time at the same rate. 

Beyond linear order, the equations of motion can be 
integrated for an Einstein-de Sitter universe, yielding for 
the growing mode (Peebles 1980), 



S2 = ^5i + Sl^a A,, 



+ -Z A,a^ A,c«/3 



(6) 



where _„ are partial space derivatives and we denote, 
A(a;) = -47r J d^X' 5{X')/\X - Xl\. 

The second order solution Eq. H, explicitly shows that, 
although $2 is of order (5f , PT is formally non-local already at 
the second order as explicitly realized by the term account- 
ing for the potential generated by the density fluctuations at 
different points (~ A). In fact, the non-local nature of grav- 
ity is an intrinsic feature of General Relativity and remains 
so in the Newtonian limit through the tidal forces when the 
appropriate limit to the equations of motion is taken as was 
recently shown by Kofman & Pogosyan (1995). 

To win simplicity in the equations of motion it is con- 
venient to work in Fourier space: 



5{k) = TT d'^s S{X) exp {ik ■ X), 



(7) 



so that one can formally write, 5{k) = ^°^ (5„(fc), being 
Sn{k) the n-th order perturbative contribution. This is ex- 
pressed as an n-dimensional integral over the kernels, Fn, 
that encode the non-linear coupling of modes due to the 
gravitational evolution, 

<5„(fc) = d''qi...d:'qiSDik-q^...-qJ 

X F„(q^...qJS^{q^)..MqJ (8) 

where 5d is the Dirac function and the kernels F„ are given 
by symmetric homogeneous functions of q^ with degree 
zero, that is, some geometrical average (see Fry 1984, Goroff 
et al. 1986, Jain & Bertschinger 1994, SF96a). In particular, 
for the second order we have: 



7 2 qiq2 \q2 qi J I \ qiq2 ) 
which reproduces Eq.M in Fourier space. 

2.3 The Monopole Approximation 

If we are in a situation where there is spherical symmetry we 
can substitute the kernels Fn in Eq. [H] by numbers Vn ■ This 



numbers are given by the monopole term in an expansion 
of the kernels in spherical harmonics as they correspond to 
their angle averages (the outcome of the angular integrals in 
Eq. l24f|). In other words, if we decompose the kernels, Fn in 
multipoles. 



Fnikl,...,kn) — Fn + / ^ -Fn i 



(10) 



only the first term F^j^", the monopole (angle average of Fn) 
which is the spherically symmetric contribution to the ker- 
nels, contributes to the Gaussian tree- level. We shall define, 
for convenience, the monopole in terms of some numbers i^n 
in the way. 



F" 



= < Fn >= Un/n\, 



(11) 



where the factor n! refiects the fact that we are lookiiig for 
some symmetric amplitudes. For instance, from Eq.|9| we 
have. 



f: 



E< Fa >= i/2/2! = 34/21. 



(12) 



On the other hand, we can write the density fiuctuation as 
formed of some generic local and non-local components. 



s{x) = s'°%x) + s'''°''{x) 



(13) 



where local means that the evolved (non-linear) density con- 
trast at a given point is only a transformation of the linear 
density contrast at the same point, what we shall call a 
local-density transformation, S''°'^{X) — £,[Si{x)]. The non- 
local component gives the contribution from density fluc- 
tuations at any other points. This is non-vanishing to all 
orde rs beyond the second in PT as commented above (see 
§2.2) since gravity induces a non-local {i.e., non-isolated) 
evolution of the density fluctuations. The equivalent to this 
local-density transformation in Fourier space would then be, 
5'°'^{k) — C[5i{k)]. In particular, we can write the n-th 
perturbative order of 5'°'^{k) in the following way. 



5l:^k) = ^ 5,{k) : 



■^Si{k) 



(14) 



which involves n convolutions of Si{k). By substituting F„ 
by its monopole contribution, Vn/nl, in Eq.Q, we immedi- 
ately see that the above Cn numbers are just those given by 
the monopole term, i.e., Cn = Vn- Transforming back to real 
space, we find Sn^iX) = Cn/n\ [5i{X)]'^ , and thus that the 
monopole approximation is given by the local contribution to 
the density fluctuation and has the form, 
00 00 

5"'^{x) = J2 ^"'(^) = ^M^)] = Y.^ MxT, (15) 

n— 1 n^l 

with ci = 1, to reproduce the linear solution, and c„ — i^n 
to recover the monopole. 



3 CUMULANTS IN PERTURBATION 

THEORY AND OTHER APPROXIMATIONS 

We will now show the connection between the tree-graphs in 
the diagrammatic approach to the one-point statistics and 
the monopole contribution to the kernels. This will allow us 
to use a simple local non-linear transformation (what we 
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shall call the local-density transformation) to derive the pre- 
dictions for the cumulants in the perturbative regime from 
the monopole contribution alone. 



3.1 Statistical Properties 

The perturbative solutions above can be used to study the 
evolved statistical properties of the density fluctuations, 
which is our final goal. Here we concentrate on the Jth- 
order moments of the fiuctuating field: mj = ( S^ Y In this 
notation the variance is defined as: 



Var{S) 



; 7712 — rrii 



(16) 



In general, it is interesting to introduce the connected mo- 
ments ^j, which carry statistical information independent 
of the lower order moments, and are formally denoted by a 
bracket with subscript c: 



^^-(^^>. 



(17) 



The connected moments are also called cumulants, reduced 
moments or irreducible moments. For a Gaussian distribu- 
tion ^,/ — for J > 2. In many cases the higher order cumu- 
lants of a particular distribution are given in terms of the 
second order one ^2- It is therefore convenient to introduce 
some more definitions. A dimensional scaling of the higher 
order moments in terms of the second order one ^2 ~ cr^ is 



given by the following ratios: 



i 



.//2 



(18) 



We will see that for Gaussian initial conditions in pertur- 
bation theory, it is more useful to introduce the hierarchical 
coefficients, 



Sj 



ij 



ii 



(19) 



as PT predicts this quantities to be time-independent quan- 
tities on large scales. These amplitudes are also called nor- 
malized one-point cumulants or reduced cumulants. We shall 
also use skewness, for Sz and kurtosis, for S4. 

The way to proceed is to substitute the PT calculations 
Eq. H and commute the spatial integrals in (S„ with the ex- 
pectation values (..■), by using the fair sample hypothesis 
(§30 Peebles 1980). Thus the calculation will reduce to av- 
erages over moments of Si with the kernels in Eq. M , which 
will, in turn, give us the statistical properties of 5 in term 
of the ones in the IC. 



3.2 Linear Theory 

If we only consider the first term in the PT series, Eq. W, and 
the growing mode, Eq. H], the cumulants of the evolved field 
will just be, (S'') = D'^ (^0) y were ( ^0 ) correspond 
to the cumulants in the IC. Consistently, the hierarchical 
ratios (see Eq.Q) will scale as, Sj = Sj{0)/D-'~^, were 
Sj(0) are the initial ratios. Note that this implies that the 
linear growth erases the initial hierarchical ratios, so that 
Sj -+ 0, as time evolves (and D —> 00). 

In terms of the dimensional scaling, see Eq. |18] , we have, 
Bj = -Bj(O), so that the linear growth preserves the initial 
values. Note that if we want to do a meaningful calculation 



of these ratios or the cumulants, in general, we might need 
to consider more terms in the perturbation series, Eq.|4|, 
depending on the statistical properties of the IC, e.g., how 
they scale with the initial variance, which is typically the 
smallness parameter in the expansion of the cumulants. 

For GIC both Bj{0) = Sj(0) = 0, and we have to 
consider higher order terms in the perturbation series to be 
able to make a non- vanishing prediction. 



3.3 The Tree-Level and Tree-Graphs 

The computation of the cumulants in PT dates back to Pee- 
bles (1980) work where the leading order contribution to the 
skewness was obtained making use of the second-order PT 
as given in Eq.|6| to give. 



S-i 



in. 



34 

7 



+ OK], 



(20) 



in agreement with the hierarchical scaling: ^j ~ ^2 • Fry 
(1984) extended Peebles analysis by making the connection 
between tree diagrams (or tree-graphs) and the perturbative 
contributions to leading order in the Gaussian case. With 
the help of this formalism he was able to obtain the lead- 
ing order contributions for the three- and four-point func- 
tions making use of the 2nd and 3rd order in PT. Later on, 
Bernardeau (1992) found the generating function of the one- 
point cumulants to leading order for GIC. Furthermore, Fry 
(1984) found that, in general, the lowest order (tree-level) 
connected part that contributes to the J-order cumulant, 
^j, is of order 2( J — 1) in 5i from terms like: 



{5lsi-')+... + {5(-'5.j-.). 



(21) 



Note that this involves the cancellation of J — 2 contribu- 
tions to the connected moment of order J. This is a property 
of the Gaussian initial conditions for which all (5(\ van- 
ish for J > 2. The above result Eq. l2lf| can be understood 
as follows. We are only interested in estimating the leading- 
order contribution (in terms of 5\) to the connected part of 
the spatial average of the product 5(k\)5{k2)...5{kn), e.g., 
that contribution made with terms that can not be factor- 
ized as products containing disjoint subsets of the labels. 
Each 5 is to be replaced with the perturbative expansion 
(see Eq. W), so that it may contain an arbitrary power in 5\. 
In the Gaussian case, the moments factorized only by pairs, 
so that the smaller contribution, with the smaller number 
of Si, connecting J points has J — 1 pairs. Each pair corre- 
sponds to two Si and therefore the leading contribution has 
indeed 2(J — 1) terms. These are the tree diagrams, corre- 
sponding to graphs with no loops (as any graph involving 
loops implies higher orders in Si). Therefore we see that for 
GIC the leading-order contribution to the cumulants (the 
so-called tree-level) is only made of tree diagrams (or tree 
graphs). That is why it is called the tree-level. We shall see 
in a companion paper (see Paper II) that the leading-order 
to the cumulants for non-Gaussian IC is not solely made of 
tree diagrams as there are additional terms which depend 
on the initial n-point correlation functions. 
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3.4 The Tree-Level and the Monopole 
Approximation 

Consider a generic term of a tree-level contribution such as 
((5„(fei)5i(fc2)...(5i(/cn-i) )• From Eq.[^ we have: 



(5„(fci)<5i(fc2)...(5i(fc„-i)) = 
X fo(fci-q,...-q„)F„(g,...qJ 



d qi...d q„ 



{5r{q,)..MqJ&i{k2)..Mk^-i)) 



(22) 



for Gaussian initial conditions. Note that the leading con- 
tribution to Sj is enough to specify all the coefficients c„ 
in the transformation, and therefore all higher order cor- 
rections. As expected, by substituting C2 by the monopole 
contribution: C2 = V2 in Eq.lld], one recovers Peebles (1980) 
calculation, Eq.|20|. 



3.6 The Local-Density Approximation and the 
Equations of Motion. 



Thus we only have to perform the spatial integral over the 
initial moments \5^ V In the case of GIC these mo- 

ments are just products of two-point functions. The only 
terms in the products of two-point functions that produce 
connected graphs have pairs connected like: 



d?qi. 



.d\^ Soik^ - q^... - qj F4q^...q^ 



X ( <5i(g,)5i(fc2) ) ( 5i(q,)5i(fc3) ) . . . ( <5i(gJ5i(fc„-i) ), 

e.g., those terms with pairs of q and k, as otherwise there 
is a part of the integral that factorizes. By isotropy, one 
sees that the only dependence in the angles between the ^'s 
comes through _Fn(<7i ...<7^). Thus, the geometrical depen- 
dence in Eq.|24[ can be integrated out and the contribution 
of the kernels Fn just becomes a constant. Thus, the (Gaus- 
sian) tree-level is given by the monopole and therefore by the 
local contribution to the density fluctuation (see §2.J). As a 
result, in this scheme, the higher-order multipoles in the 
kernels are only expected to yield a non- vanishing contribu- 
tion in the next-to-leading orders (loops) in the cumulants, 
while the monopole contributes to all perturbative orders. 
In what follows we shall deal with the local contribution to 
the density ffuctuation alone, Sioc (monopole in the kernels) 
and work out its 1-point statistics. 



We shall stress that the local-density transformation (see 
Eq.llql) is just a particular case of what is usually under- 
stood as local. In general, it means that the value of the 
evolved field at a given point is a transformation of the IC 
at the same point, say, the initial value of the field, its deriva- 
tives and any other fields {e.g., the velocity field) contribut- 
(23) ing at that very point. Although gravity is non-local in this 
restricted sense, the above results show that the statistical 
average of the non-local integrals involved in the reduced 
moments to leading order are exactly given by a local- density 
transformation. 

However, to estimate the leading order (tree-level) it is 
not necessary to calculate the full kernels F„, as we only 
need the c„ numbers, the monopole (spherically symmetric) 
contribution. Therefore, the values of Cn (and thus the ap- 
propriate local-density transformation) can be determined by 
finding the spherically symmetric solution to the equations of 
motion. We shall see below that for the gravitational evolu- 
tion of the density field in Lagrangian space, this transforma- 
tion is given by the spherical collapse (SC) model, whereas 
for the ZA to the dynamics, it is given by what we shall call 
the Spherically Symmetric ZA (SSZA). 



3.5 Estimation of the Cumulants 

We can now easily estimate all the 1-point statistical prop- 
erties in the monopole approximation to PT. This can be 
done by using the generating function method: 



^^ = {S'').= 



d-' ln[Z] 



dtJ 



(24) 



where Z — (e^ ^ is the moment generating functional of the 
one-point probability distribution functional, and the field 
5 = Sloe is given by Eq.[tl5|. Furthermore, ln[Z] is the cu- 
mulant generating functional so that Eq.|24| yields a simple 
derivation of the cumulants in terms of the moments. A gen- 
eral discussion in terms of the p-point correlation functions 
and its computation using diagrammatic techniques (simi- 
lar to those in Quantum Field Theory) is given in detail in 
SF96a. 

The resulting expressions can be found in Fry & 
Gaztafiaga (1993, FG93 hereafter), who consider a generic 
local transformation between density fields and find, to lead- 
ing terms in ai : 

Ss = 3C2 + 0{af) 

Si = 4C3 + 12ci -I- 0(af) 

Ss = 5c4 + 60C3C2 + 60ci 4- 0{af) 



3.7 Other Approximations 

In one dimension, the probability distribution function of the 
matter field may be integrated in Lagrangian coordinates 
making use of the Zel'dovich Approximation (ZA) which 
recovers the exact dynamics. For the three-dimensional case 
however, the ZA fails to describe the exact picture since the 
displacement field (that only depends on the IC) no longer 
factorizes in the actual solution for the growing mode as the 
ZA assumes (Bernardeau & Kofman 1995). 

In the context of Lagrangian space and inspired by the 
successful results of the Zel'dovich Approximation (hereafter 
ZA), there have appeared in the last years a number of pa- 
pers presenting some local approximations to the exact dy- 
namics that successfully describe the evolution of the den- 
sity contrast in its first non-linear stages (see Bertschinger 
and Jain 1994, Hui and Bertschinger 1996). Most recently, 
Protogeros and Scherrer (1997, hereafter PS97), following 
the formalism of the ZA, made use of a closed analytic ex- 
pression, 5 — (1 ~ Si/a)~°' — 1, to derive the cumulants 
following the formal analogy between the density contrast 
and the vertex generating function at tree level (see B92), 
for different values of a, which they take to be different ap- 
proximations to the exact PT. For the particular case when 
Q = 3 (what they call the "spherical approximation") their 
approximation happens to be equal to the spherically sym- 
metric solution to the dynamics (in 3 dimensions) for the 
(25) ZA (what we call SSZA in Appendix O). They also consider 
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the case a = 3/2 as another approximation they call the 
"exact approximation" (which, in fact, only reproduces the 
exact tree- level when f2 = 0). In general, none of these local 
approximations recovers the exact tree-level amplitudes as 
predicted by PT (which we will show later that is given by 
the SC model). 

On the other hand, a number of non-local approxima- 
tions to the non-linear dynamics of the cosmological fields 
have been designed as to simplify the equations of motion 
by linearising any of the fields involved. The Linear (or 
Frozen) Potential Approximation (LPA) which assumes that 
the gravitational potential generated by the density pertur- 
bations remains linear throughout the evolution (see Brain- 
erd et al. 1993, Bagla and Padmanabhan 1994), the Frozen 
Flow Approximation (FFA) that takes the velocity field to 
be frozen to its initial shape (see Matarrese et al. 1992), 
along with the ZA, which deals with the approximation that 
particles move on straight lines in the comoving picture (see 
Zel'dovich 1970), are the best known among them. A general 
treatment in terms of the vertex generating function can be 
found in Munshi et al. (1995). For a review of all these ap- 
proximations see Bernardeau et al. (1994). If we write the 
second perturbative order as given by Eq.lq], in the generic 
form: 



5i +A5i + B5i,c/^,c + CA,c«/3A,c/3, 



(26) 



one may derive the skewness at tree-level in the parametric 
fashion, 



Ss = 9^-3B -f 5C, 



(27) 



from which we can characterize the different above cited ap- 
proximations depending on the weight each of them gives to 
these three different contributions: while PT gives A = 5/7, 
B = 1, C = 2/7, the ZA yields 1/2, 1, 1/2, whereas the FFA 
gives 1/2, 1/2, 0, and the LPA gets 1/2, 7/10, 1/5, as the 
corresponding values for A, B and C respectively, in each 
approximation. By just replacing the numbers in S3 one sees 
that none of these non-linear approximations is able to re- 
produce the exact tree-level amplitude given by PT for the 
skewness. This illustrates that, despite successfully describ- 
ing the first non-linear stages of the gravitational collapse 
of a density fluctuation, none of the above given approxima- 
tions yields to accurate enough predictions concerning the 
one-point cumulants in the perturbative regime. 

3.8 Higher Perturbative Orders 

Whenever 5 < 1, corrective terms beyond leading order 
(tree-level) become quantitatively important and must be 
taken into account. A systematic approach for deriving the 
next-to-leading orders in PT was introduced by Goroff et al. 
(1986) by representing them in terms of the loop diagrams. 
Recently, the first results for the loop corrections in PT have 
been calculated (see SF96a, SF96b and S97). 

The angle averaged (spherically symmetric) picture de- 
scribed above can be used to estimate the monopole contri- 
bution to an arbitrary higher-order in PT. This is easy to 
do as one just has to take more terms in the series given 
by Eq. |1^ to calculate the next orders in the moments. As 
expected, even in the Gaussian case, the higher-order terms 
coming from the monopole contribution differ from the exact 
PT estimates. Bernardeau (1994c, B94c hereafter) already 



realised that fact when tracing the non-linear evolution of 
a volume occupied by particles that were in a given initial 
density perturbation. If one assumes that the matter content 
of a fluctuation is conserved during its collapse, then tracing 
the evolution of the volume is equivalent to following that 
of the matter density. In his work, the density fluctuation 
follows the SC model in what he called the rare event limit, 
I.e., when Si/ai — > 00, thus it is only exact as a mean picture 
of the actual collapse, with no scatter, i.e. oi = 0. 

We have seen above that, in Fourier space, the depar- 
ture of the SG picture from the exact PT is because the 
exact integrations in this case, e.g., such as that in Eq.|22| 
with more 5i factors, involve loops and the geometrical de- 
pendence is not trivial anymore (the dipole, quadrupole... 
contributions no longer cancel); The kernel Fn integration 
will give a different number depending on the loop configu- 
ration (see SF96a), contrary to the monopole approximation 
(were F„ is always replaced by the same amplitude Cn/n\). 
However, one can naively expect the monopole contribution 
to dominate by symmetry because asymmetric contributions 
(arising from tidal terms) tend to cancel when averaged. In 
this paper we shall compute the monopole contribution to 
PT and will compare it to the exact analytic results, when 
available, or N-body simulations, to see how accurate this 
approximation is. 



4 THE SPHERICAL COLLAPSE MODEL 

4.1 The Shearless Approximation in Lagrangian 
Space 

We next want to derive the evolution of a density fluctua- 
tion within the spherical model of collapse. In order to do 
that, we turn to the Lagrangian space which is the natural 
framework for describing the motion of a fluid element. The 
SG dynamics is fully described in terms of one parameter, 
the initial size of the spherical fluctuation. In this sense, the 
SG model is just a particular case within the family of local- 
density approximations discussed above. When we write the 
equations in Lagrangian space, the natural variable is the 
density contrast and the only parameter in terms of which 
the solution is given is the linear density contrast. 

We first must define the conformal time r which is the 
comoving time parameter to the motion of the mass and 
we shall denote with a dot the associated time derivatives 
(conformal time derivatives), defined as. 



d d „ 

dr or 






(28) 



a(t) being the scale factor. The generic equation describing 
the evolution of a spherically symmetric perturbation in an 
expanding universe is given by. 






GMjR) 

i?2 



= -A-KGpR 



(29) 



since the matter contained in a spherical perturbation of ra- 
dius R is, M{R) — A-RpE? /i. Therefore, the spherical den- 
sity perturbation 5 = (R/a{t)Ro)~ — 1, is described by the 
following equation of motion 



5 + H{r)5 



3 (1 + 5) 



47rG(p-p)a(r)2(l + <5). (30) 
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where the dot denotes a conformal time derivative, and a^r) 
is the scale factor in terms of the conformal time from 
which we define the conformal Hubble parameter 7i(r) = 
d log a(r)/dr. Introducing the comoving time derivative, 
d/dt = ■, the last equation translates into, 



S +2H{t)5 - - 



3 (1 + 5) 



= 47rGp5(l + 5), 



(31) 



where H{t) = d\oga{t)/dt. 

On the other hand, the Newtonian Fluid equations with 
zero pressure, which are the appropriate description for the 
sub-horizon modes in a perturbed FRW universe, in its mat- 
ter dominated regime (the relevant for describing formation 
of structures), are usually given in Euler space (see e.g., §9 
in Peebles 1980). If we now turn to Lagrangian coordinates 
and we re-express the equations of motion in terms of the 
derivatives of the conformal time r, the continuity equation 
reads, 



(1 + 5)6* =0, e= V -v. 



(32) 



We then combine this last expression with the Raychaudhuri 
equation. 



+ n{r)e + -( 



+ cr'V, 



2lj^ = -A-KG-pSa"^ 



(33) 



where to = ■^uj'-'lo, 



ij, and the expansion 9, vorticity ujij, and 
shear aij , are given by the trace, traceless antisymmetric and 
symmetric parts respectively, of the velocity divergence. 



ViVj 



9Sii + a-ii + uji 



(34) 



and get a second-order differential equation for the density 
contrast, 



5 + n{r)5 



3 (1 + 5) 



2lu'^ + A-jvG-pSa 



(35) 



For an initially irrotational fluid (the expansion preserving 
its irrotational character in the linear regime): oj = 0. Mak- 
ing the further assumption that there is no shear, Eq.|35| 
leads to the equation we obtained for the SC model (see 
Eq. |3C[| above). In other words, the SC approximation is 
the actual dynamics when tidal effects are neglected. As one 
would expect, this yields a local evolution, in the restricted 
sense that the evolved field at a point is just given by a local 
(non-linear) transformation of the initial field at the same 
point. Throughout this paper we shall drop the shear term 
and work out the solution to this dynamics in the perturba- 
tive regime. 

The exact (non-perturbative) solution for the SC of the 
density contrast in an Einstein-deSitter universe admits a 
well-known parametric representation. 



^(0)4-!^^ 



2(1 



SiW = ^ 



- (0 - sin 4>) 



2/3 



for Si > 0, linear overdensity, and 
9 (sinh0-0)2 



SW 



2 (cosh </>- 1)3 



SiW 



-(sinh( 



2/3 



for 5i < 0, linear underdensity (see Peebles 1980), where the 



parameter (j) is just a parameterization of the time coordi- 
nate. 

If we are only interested in the perturbative regime 
{5i — > 0), which is the relevant one for the description of 
structure formation on large scales, the above solution can 
be also expressed directly in terms of the initial density con- 
trast, which plays the role of the initial size of the spherical 
fluctuation in Eq.pl|. This way, the evolved density con- 
trast in the perturbative regime is given by a local- density 
transformation of the linear density fluctuation. 



/w-E^[^'] 



(36) 



Notice that all the dynamical information in the SC model 
is encoded in the Vn coefflcients of this local-density trans- 
formation, Eq.fsd]. In an Einstein-de Sitter universe (fi — 
1, A = FRW universe), we can introduce the above power 
series expansion in Eq. IM and determine the u„ coefflcients 
one by one. The flrst ones turn out to be. 



2^2 



34 
21 
446440 



1.62; 



1^3 = 



682 
189 



U4, = 



10.22; 



-3.61 
8546480 



r^s = 



35.13 



(37) 



43659 ' 243243 

and so on. 

Remember that this evolution is in Lagrangian space, 
5 = S{q). We would like to relate the above results, ob- 
tained in Lagrangian coordinates to the corresponding fluc- 
tuation in Eulerian coordinates. In Lagrangian space a fluid 
element of a given mass is labeled by its initial position q 
(or Lagrangian coordinate), whereas Eulerian space uses the 
density p{x) at the final coordinate X — X(q). Notice that 
mass conservation requires that the volume elements be re- 
lated like: d^q = (1 + 5)d^X. Thus density probabilities in 
Lagrangian (_l) and Eulerian (_e) space should also be re- 
lated by the same factor, AlS = (1 + S)AeS (see Kofman 
et al. 1994). This provides a simple way to translate density 
moments in Eulerian space, ((1 + 5)'^ ) with the ones in 

Lagrangian space ( (1 + 5)'^ ) , 



((1 + 5)'')l= {l + syP[5]AL5 = 



(1 + Sy+^P[5] AeS = ((1 + 5y+')E. 



(38) 



In particular, as pointed out by Bernardeau (B94b) (and also 
PS97), the conservation of Eulerian volume yields, (1)_b = 
1 — (( 1 + 5)~^)l, which requires a normalization for 5 in 



l + 5 = (l + /)((l + /)-i)^. 

It is interesting to note from Eq.bsl that: 

{S')l^{S')e + {S-'+')e. 



(39) 



(40) 



Thus, to leading order, there is no difference between Eule- 
rian and Lagrangian moments. In general, the leading con- 
tribution to the cumulants is enough to specify all the co- 
efficients c„ (see e.g., Eq.[G5[), of the local transformation 
provided by Eq.|l^. This means that the local-density trans- 
formations that produce the tree-level in PT are identical in 
Euler and Lagrangian space. A similar result was noted by 
Protogeros & Scherrer (PS97), in the context of hierarchical 
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distributions, but this argument is more general as applies 
to any non-Gaussian distributions where: {5'^^^)e < {S'')e, 
which holds even for the strongly non-Gaussian dimensional 
models, Eq.|tL8|, that we shall investigate later (see Paper 
II). These models for the IC naturally appear as solutions 
for the topological defects models for structure formation. 



4.2 The SC as a Solution to the Cumulants 

Our goal here is to relate the cumulants of the initial dis- 
tribution with the ones resulting from the full non-linear 
evolution of the field. To estimate the evolved cumulants 
under the SC model, we perform the following prescription: 

• (a) we start with the initial cumulants, in the limit of 
small fluctuations, 5 -^ . In this limit, the cumulants are 
equal in Lagrangian and Eulerian space, as, d^q ^t d^X. 

• (b) we then use the local transformation Eq. [M, with 
its proper normalization, Eq.[p9|, to relate the initial and 
final cumulants in Lagrangian space. As argued above, this 
is equivalent to include the contribution from the monopole 
alone. This can be done as described in section §B.5 using 
Eq. |24] to obtain expressions such as Eq. psl] by keeping the 
relevant terms. 

• (c) finally we use Eq.lSSl to relate the Lagrangian and 
Eulerian moments, which can be rewritten in a more com- 
pact way: 



{5')e = {5'-')l-{5''^)l + {5'''')l 



+ (5)i, (41) 



from which we get the cumulants in Euler (real) space. 

Note that the cumul ants obtained from (b)-|-(c), which are 
Appendix M (see Eq. A2|), are not the same as the ones ob- 
tained using the local-density relation in Euler space directly 
(see Appendix H). As explicitly seen from Eq.[J40|, they are 
only the same to leading order. We choose this approach 
because the Lagrangian space, comoving with the fluid ele- 
ment, is the natural coordinate system for a local description 
of its evolution. Note that a local-density description in La- 
grangian space will in general mean a non-local one in Euler 
space. 



4.3 The SC Model and PT 

We have pointed out in the previous sections that the 
monopole contribution (which is the exact result for tree- 
graphs) to the cumulants in PT is given by a local-density 
transformation Eq.||l5|, whose coefficients, c„, are to be de- 
termined by the kernels, Fn (see Eq. S), under the relevant 
dynamics. These arguments are valid for any dynamics and 
apply either to Euler or Lagrangian space. They are true for 
any leading order calculation. As argued in section ^3.6| , to 
estimate this contribution to the cumulants it is not neces- 
sary to calculate the full kernels F„, as we only need the 
numbers c„ . Given the equations for the evolution of a field, 
one can determine Eq.[[15| and therefore c„ by just requiring 
the solutions to be spherically symmetric. 

We have shown above that for gravity, the spherically 
symmetric solution to the evolution of density perturbations 
is given by Eq.[^, e^., the SC model, whose solution is 
well-known (see Eq.fSq]). Thus, £, — f which yields the 



monopole contribution, c„ — i^„, without need of estimat- 
ing the kernels F„ or any integral. Note, in particular, that 
U2 = 34/42 from both approaches (see Eqs.jl2|,||^). 

This provides with a simpler derivation and interpreta- 
tion of the results presented by B92, who found the values i^n 
that give the leading-order contribution to PT for Gaussian 
initial conditions. In the language of B92, the vertex gener- 
ating function, G{—t) = f{Si) (see also PS97). Our deriva- 
tion explains therefore why the vertex generating function 
Q(—t) follows the SC model, which lacked a satisfactory ex- 
planation in the context of B92. Besides its simplicity, in 
our framework one has the added advantage of being able 
to use the local-density relation to estimate the higher-order 
corrections for both Gaussian and non-Gaussian initial con- 
ditions (see Paper II). Note nevertheless that there is an 
important difference in practice with the work of B92. His 
vertex generating function Q{—t), corresponds to cumulants 
in Euler space, while our local-density relation f{St), applies 
to Lagrangian space. To leading order, both give identical 
results for Gaussian IC, but they yield different results in 
general for Non-Gaussian IC or for next-to-leading terms 
for Gaussian IC. 

We want to stress that the fact the SC model deter- 
mines the tree-level amplitudes is neither a singular prop- 
erty of the Gaussian nature of the IC nor an argument that 
has to do with the flat geometry of space. It is a general 
property that follows straightforwardly from the local na- 
ture of the monopole contribution to PT (see section §3.6): 
the tree-graphs (and the monopole contribution) in PT are 
exactly given by the SC model regardless of either the sta- 
tistical nature of the initial conditions or the geometry of 
the universe. We shall illustrate the above statement in the 
framework of B92 calculations by showing, in an accompa- 
nying paper (Fosalba & Gaztanaga 1998, Paper III in this 
issue) , that the equations of motion that govern the leading 
order amplitudes (for Gaussian initial conditions) in a uni- 
verse with arbitrary density parameter fi, are those of the 
SC dynamics. 

4.4 Smoothing Effects 

So far we have worked out the statistics of the unsmoothed 
fields. There remains to be seen whether the local-density 
approximation may be extended to the statistics of the spa- 
tially smoothed fields, which is essential if we are to compare 
with N-body simulations and observations. For that pur- 
pose, we have to introduce the fluctuating field integrated 
over a finite volume. This volume is fixed by the size of the 
window function that acts as a filter on the unsmoothed 
field. 

Our goal is therefore to relate the smoothed cumulants 
in the evolved distribution with the smoothed cumulants in 
the initial one. The latter are inputs to our calculations 
and should tell us, e.g., how the smoothed rms fiuctuation 
changes as a function of the smoothing radius: a — a{R). 

We will focus here in the top-hat filter, defined as. 



WTHiX,R) 



1 if X <R, 



being zero otherwise, where R is the smoothing radius. Note 
that this filter is spherically symmetric so that the trans- 
formation that gets from the unsmoothed to the smoothed 
evolved flelds preserves the spherical symmetry. 
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In the local-density picture of the spherical collapse, 
each fluctuation is isolated from the others and evolves ac- 
cording to the value of the initial amplitude alone. The 
statistics of the evolved field is just induced by the statistics 
of the initial one, which could in general be non-Gaussian. In 
order to estimate the statistical properties (one-point cumu- 
lants) we can therefore picture the initial spatial distribution 
as just made of a single spherical fluctuation. Different real- 
izations of this fluctuation could have different amplitudes 
or shapes in a proportion given so as to match the statistics 
of the initial conditions (see §61 Peebles 1980 for examples). 
In our case this is given in terms of the unsmoothed distri- 
bution, which corresponds to a top-hat smoothing of radius 
Ro -^0. 

The spatial smoothing only acts on the fluctuation by 
changing its amplitude. We want, by construction, the new 
amplitude to be a function of the smoothing radius, R. Thus, 
for the initial conditions Si, the smoothed field is: 

m) - ^§^A (42) 

cr;(ito) 

as smoothing does not change the statistics directly (at least 
for a top-hat). A similar argument applies to the evolved 
field, which under the SC model is just a local transforma- 
tion of the initial fluctuations, which only changes the local 
amplitude: 



S^^S 



(43) 



where S is the smoothed fluctuation and 5 the unsmoothed 
one. The statistics of 5 will therefore be given by the un- 
smoothed statistics, S, which are in turn given by the local- 
density relation Eq.Q: 5 ~ 5 = f[Si]. Using Eq.H we 
express this as a function of the linear smoothed fluctuation 
Si: 



5[5i]-m-^] 



Ol 



(44) 



where ~ just means that 5 needs to be normalized, as ex- 
plained before (see Eq.|39|). Note that the input rms fluc- 
tuation is a function of the smoothing radius, R, which also 
fixes the amplitude of the final smoothed fluctuation. As the 
smoothing radius R is, by construction, larger than the un- 
smoothed one i? >> _Ro — * 0, we have for the smoothed 
fluctuation: p — ^ ~ mo R~'^, as there is no other smooth- 
ing than the sharp cut-off, (note that this may not be the 
case for other filters, such as the Gaussian smoothing, which 
changes the initial shape of the fluctuation and therefore the 
final mass). Thus the smoothed amplitude, a — o-[Ki will 
be a function of 5 through R. This together with Eq.|U4| can 
be used to obtained 5 in a recursive way. 

In the case of a power- law power spectrum P{k) ^ k" , 
the smoothed variance is also a power-law ai ~ i?^", where 
■y = -{n + 3). We then have, en = a; (1 -I- Sy'^. Note that 
7 = reproduces the unsmoothed result. Moreover, using 
Eq.@, we find: 

m-mi^ + sy^'] (45) 

up to a normalization factor given by Eq. pa] . Note that this 
final result as well as the general expression Eq.|^, agrees 
with B94a arguments, based on the vertex generating func- 
tion. However, Eqs.p4|,|4a] do not limit themselves to Gaus- 
sian initial conditions or the leading-order term. Here again. 



the vertex generating function Q{~t), corresponds to cumu- 
lants in Euler space, while our local-density relation, f{5i), 
applies to Lagrangian space. To leading order, they both give 
identical results with Gaussian IC, but they yield different 
results, in general, for Non-Gaussian IC or for higher-order 
terms with Gaussian IC. 

We will further write the smoothed density contrast 
(over a smoothing scale R) in the following way: 



^ = /w = Ef[^'] 



(46) 



where the Vk are a generalization of the unsmoothed coeffi- 
cients. By Taylor-expanding Eq. [ttsl one obtains. 



1^2 — V2 + 



7 



1 2 

V'i = -(-27 + 7 +67!/2+4!/3) 

W = — (367-367^ + 87^ -1087!/2 + 72 7^ !/2 + 

+ 5471^1 + 7271^3 + 271^4), 



(47) 



the above results 
The results are 



and so on. 

For an arbitrary power spectrum P 
can be trivi ally generalized usingEq.| 
given in Eq. ]A4| of the Appendix ^. 

In Appendix O we present an alternative derivation of 
the (top-hat) smoothing effects within the SC approxima- 
tion, following the kernels in Fourier space. There, one can 
explicitly see how in the SC model, smoothing is a trivial 
operation which only changes the Vn coefficients but not the 
local nature of the transformation from the linear to the 
non-linear density field. 



5 COMPARISON OF THE SC MODEL WITH 
EXACT PT AND N-BODY SIMULATIONS 

In this section we start from the exact solution to the SC 
dynamics in an Einstein-de Sitter universe assuming the rms 
fiuctuation ai < 1. In that perturbative limit, we carry out 
the calculation of the connected moments for the density. 
Results for the velocity field are provided in Paper HI. 

We compare the SC predictions with those derived from 
the exact PT in the context of the diagrammatic approach. 
Although we know that the SC dynamics exactly reproduces 
the leading-order contribution to the cumulants, higher or- 
der effects will in general be different. We have seen that the 
difference is due to neglecting the shear and can therefore 
be attributed to non-local effects, which we will also call 
tidal effects. Note that in our context local refers only to the 
density. 

As commented above, the results are derived in gen- 
eral for the smoothed fields for a top-hat window function, 
the unsmoothed fields being recovered just as the particular 
case 7 = 0. To simplify the expressions, we will assume that 
higher order derivatives of the variance (see Appendix Rl 
Eq. A4| with 7^ ~ for p > 2) can be neglected. This is a 



good approximation for slowly varying P{k), like CDM (see 
Scoccimarro 1998) or the APM model, but it is straight- 
forward to take those higher order derivatives into account 
anyway. We do take into account 72, which for loop correc- 
tions could contribute up to 10% on small scales. 
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We will also compare the SC results with different N- 
body simulations, with parameters given in Table [U Also 
given is the reference where the details for a particular run 
can be found. For each simulation, the cumulants ^ j are esti- 
mated from counts in spherical cells, as described in BGE95, 
where more details about the estimation are given. For the 
later outputs the correlations are corrected from the Poisson 
shot-noise {e.g., BGE95). We also apply the corrections due 
to the ZA transient (see Baugh et al. 1995) as analitycally 
derived by Scoccimarro (1998). 



4.2 



Following the steps described in section 
derive now the smoothed one-point cumulants to an arbi- 
trary perturbative order for the SC dynamics. In this case 
we start from Gaussian initial conditions and use the local 
transformation for the local smoothed density (see Eq.M]) 
with Cfc = Vk, as given in Eq.|U7|, to find the leading-order 
and higher-order corrections for the variance and the hi- 
erarchical amplitudes. In order to handle the perturbative 
expansions in the cumulants for both GIG or NGIC in a gen- 
eral framework, we shall introduce the following notation. 



={n=Y: 



S2,i CTi 



(48) 



where S2,i = 1 and the subscript i in the coefficients la- 
bel the perturbative order. For GIG the odd terms in the 
perturbative expansion vanish and the get for the variance. 



2 2 

0"G = ^l 



4 , 6 , 

S2,4 O-i + S2,6 CTi + 



(49) 



Note that our notation for GIG is equivalent to that of SF96a 
provided one identifies S2,2i-i-2 = s^'' with i = 1, 2, • • •. 

On the other hand, for the hierarchical amplitudes we 
keep the above introduced notation with the added labeling 
of the order of the moment J, that defines the Sj coefficients. 



Sj = 



?.. 



i^'). 



-d- 



(n 



T=E^^>' 



<^i 



(50) 



which for GIG has non-vanishing contributions from the 
even terms alone. 



Sj — Sj,0 + Sj,2 CTi + S'j,4 at + 



(51) 



with Sj9i 



Sf (with 



As mentioned before. 



= 0, 1, 2, ■ ■ ■) in SF96a notation. 

it is usual to denote tree-level to 

:.(0) 



the leading order contributions, S j . Although this is true 
for GIG, in general the leading order contributions do not 
correspond to the tree-level in the diagrammatic approach 
(see Paper II). However, for the sake of clarity, we shall 



j{0) 



.(0) 



maintain the notation S''J'' (or T} ' for the velocity field, see 
Paper III) to denote the leading-order contribution. 

With this notation the results for GIG (in te rms of the 
SG coefficients Uk) are given in Appendix N, Eq.|A2 , with 
Cfc = Vk- In terms of the smoothing index 7 - 
we have, for a power-law power spectrum P{k) 
top-hat window: 

1909 143 11 2 

■7+^7 



-(n-f 3), 
fc", and a 



S2,4 



S2,6 ~ 



1323 ^ 126 

344439415 

+ 



36 
21651395 



+ 



107270163 
1651 3 



1134 



7 + 



3667356 
127 4 
648'^ 



7 + 



408721 2 
95256 '^ 



+ 



Ssfi 



5^ = 4^ + 7 



1026488 12862 407 2 10 3 
^'^ ^ 101871 "^ 1323 '^"'^ 126^ ^27'^ 

251978977148 71492200235 138567091 2 



-f 



5256237987 750891141 
79295 3 , 2891 4 , 1841 5 



1833678 



■7 +■ 



ft.o 



2646 ' 486 
(0) _ 60712 , 62 



7 + 



3888 






22336534498 42649448 3571621 2 



+ 



83432349 130977 
35047 3 1549 4 



23814 



■7 + 



■7 



1134 ' 648 
126152927186426522 69638296109567 
'*''' ^ 61923739724847 "^ 15768713961 "^ 
8977285860007 2 7018518515 3 

■7 + ogg-TOcr. 7 



+ 



2252673423 ' 3667356 
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916839 
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■7 +■ 
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■7 + 
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27810783 
138730 3 1210 4 

+ ,„^ 7 + -t:=-^ 



305613 



■7+' 



3969 



189 



27 



93347762463213320 2034255356621746 



+ 



421249930101 
211757079765188 2 



750891141 



7 + 



5256237987 
301575001360 3 
2750517 ''' 



69888305 4 1582214 5 14591 g 
+ — TTTTTT^ — 7 H ;r;^ — 7 +7 



2916 



567 



108 



(52) 



In table Q we display a summary of results for the SG 
model for different values of the spectral index. Note, to start 
with, that for large n or more negative 7 the coefficients de- 
crease with the perturbative order (the subscript after the 
coma), indicating possible convergence of the perturbative 
series. For 7 = {i.e., n — —3) the coefficients of each or- 
der in the expansion increases quickly with the order, which 
might indicate that the PT expansion does not converge for 
the unsmoothed field, at least in the SC approximation. 

For Sj, there is a suppression of non-linearities with the 
effects of smoothing (as n increases from —3 to 0), which is 
also found in numerical simulations (see e.g., Lokas et al. 
1995)). In particular, vanishing non-linearities are found for 



5.1 Comparison with Exact PT 

The tree-level results in the SC model are of course identical 
to the ones estimated in the exact PT by Juskiewicz et al. 
(1993) and B94b, as argued in §2.2 . For the unsmoothed 
fields, analytic results including the first corrective term to 
the tree-level, were first derived by SF96a through loop cal- 
culations in the diagrammatic approach to the exact PT. 
For the variance, skewness, and the kurtosis (the last one 
only at tree-level) they obtain. 



a^ ^af + 1.82 af + O (af) 
S3 ^ 4.86 + 9.80 af + O (af) 
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Table 1. Simulation parameters 



run names 



P(fc) 



n-A 



number mesh ^6oi 

of particles (h^^Mpc) 



Reference 



SCDM (a)-{j) 


r = 0.5 CDM 


n = 


= 1 A = 


126^ 


1283 


378 


LCDM (a)-(j) 


r = 0.2 CDM 


n = 


= 0.2 A = 0.8 


1263 


1283 


378 


APMPKl (a)-(e) 


APM 


n = 


= 1 A = 


1263 


1283 


400 


APMPK2 (a)-(b) 


APM 


Vl = 


= 1 A = 


2003 


1283 


600 



Gaztafiaga & Baugh 1995 
Gaztaiiaga & Baugh 1995 
Gaztafiaga & Baugh 1997 
Baugh & Gaztaiiaga 1996 



sc 


Unsmoothed 




Smoothed 






7 = 


7 = -l 


7 = -2 


7 = -3 




n = -3 


n = -2 


n = -1 


n = 


S2,4 


1.44 


0.61 


0.40 


0.79 


S2,6 


3.21 


0.34 


0.05 


0.68 


•S'3,0 


4.86 


3.86 


2.86 


1.86 


•S'3,2 


10.08 


3.21 


0.59 


-0.02 


■S'3,4 


47.94 


3.80 


0.07 


0.06 


S4,0 


45.89 


27.56 


13.89 


4.89 


S4,2 


267.72 


63.56 


7.39 


-0.16 


54,4 


2037.2 


138.43 


1.99 


0.31 


•S'5,0 


656.31 


292.35 


96.50 


16.52 


S5,2 


7242.2 


1263.97 


91.85 


-1.16 


•S'5,4 


80903.0 


4363.92 


42.89 


1.53 


Sefi 


12653.49 


4141.58 


876.62 


67.81 


Se,2 


221597.1 


28256.19 


1274.38 


-8.04 


S'6,4 


3405857.8 


140641.3 


906.74 


7.88 



Table 2. Values for the higher-order perturbative contributions 
in the SC model for the unsmoothed (n = —3) and smoothed 
{n = —2, —1,0) density fields, for a top-hat filter an a power-law 
power spectrum. 



54 ^ 45.89 + O (o-f ) 



(53) 



for the average values in the range 2 > n > — 2 (for a 
power-law power spectrum, P{K) ~ k"), with a 3% vari- 
ation within that range due to non-local effects. The above 
results are to be compared to those from the SC in the per- 
turbative regime (truncated at the same order), 

a^^af + 1.44 af + O (af) 

S3 ^ 4.86 + 10.08 0-? + O {af) 

S4 ~ 45.89 + 267.72 af + O (erf) . (54) 

Comparing with the exact PT result Eq.|5d], we can see 
that tidal (non-local) effects only amount to a 3% in the 



corrective term for Ss and up to 20% in the first corrective 
term for the variance. The 54 correction (~ 268 CT; ) must 
be taken as an accurate prediction for 54 (within the few 
per cent effect expected from the net tidal contribution for 
the Sj ratios) since there are no analytic results available to 
compare with. We therefore see from the above comparison 
that the shearless (SC) approximation leads to very accurate 
predictions for the hierarchical amplitudes, Sj, while giving 
worse estimates for the cumulants, ^j. 

We stress the importance of applying the SC approxi- 
mation in Lagrangian space, where is the exact spherically 
symmetric solutions of the field equations. The SC model 
is described by a transformation that only depends on the 
value of the linear field at the same point (what we call 
a local-density transformation). When going back to Euler 
space the density fluctuation (defined at a point) is nor- 
malized with a factor (see Eq.flSSl]) which is a function of the 
(non-linear) variance. Since the variance is a volume average 
of the two-point correlation function, this factor yields some 
non-local contribution to the cumulants (in Euler space). 
This non-local contribution is missed when introducing the 
SC model in Euler space directly, thus is not surprising that 
the predictions for the cumulants in the SC approximation 
in Euler space are a poor estimation of those in exact PT, as 
the latter are dominated by the non-local (tidal) effects (see 
Table ^ in Appendix |b|). The contribution to the cumulants 
of this non-local term is typically negative so the predictions 
from the SC model in Euler space generically overestimate 
those in Lagrangian space. However, the domination of non- 
local effects in the cumulants for the SC model in Euler 
space is partially canceled when computing the hierarchical 
ratios Sj, similarly to what is found for the SC model in 
Lagrangian space. We shall see below that smoothing effects 
do not alter substantially this interpretation (at least for a 
top-hat window). 

For the smoothed fields, there are few analytic results 
which concentrate on the variance (and its Fourier trans- 
form, the power-spectrum), the skewness, and the bispec- 
trum. For the present analysis we focus on the results for 
the variance and the skewness (see SF96b and Scoccimarro 
1997, hereafter S97, respectively) which, in the diagram- 
matic approach to the exact PT to one-loop order, for a 
top-hat smoothing and n = —2, give 



53 



3.86 



0.88 erf 
3.18 erf 



(-0 



O"! 



(55) 
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Figure 1. Departures from the tree- level contributions as the lin- 
ear rms fluctuation grows, for the variance, skewness and kurtosis 
of the density field as predicted by the SC model up to the 2nd 
non-vanishing perturbative contribution (one-loop) for different 
values of the spectral index: the dotted line shows the n = —3 
(unsmoothed) case, and the short-dashed (n = —2), long-dashed 
(n = —1), dot short-dashed (n = 0), dot long-dashed (ra = 1) 
depict the behavior for the smoothed density field. The solid line 
shows the tree-level values as a reference. The corrective term has 
a minimum contribution to the variance for n si — 1, while the 
hierarchical amplitudes show a vanishing one-loop contribution 
for n ^ 0. 



that we compare to the SC predictions, 

a^^af + 0.61<7f + 0{af) 
Sa ^3.86 + 3.21 erf + O (af) 



Figure 2. Same as Fig hi when the 3rd perturbative contribution 
(two-loop) is included in the computation of the cumulants. 



(56) 



which means a 30 % and 1 % contribution from the tidal 
forces for the corrections in the variance and the skewness, 
respectively. These can be obtained by just subtracting the 
shearless (SC) contribution to the exact computation car- 
ried out by means of the loop calculations. Thus, we see 
that in line with the unsmoothed predictions, the shearless 
contribution completely dominates the skewness S3 , at least 
at the one loop order. 

In Fig hi we display the departures from the tree-level 
contribution for the evolved variance when the one-loop cor- 
rection is included, skewness and kurtosis as the linear vari- 
ance approaches unity, where the perturbative regime is ex- 
pected to break down. Fig Q shows the same as Fig hi when 
the 3rd order (3rd non- vanishing contribution) is included in 
the computation of the cumulants: the qualitative features 
are preserved when this higher-order contribution is taken 
into account with respect to those inferred from the 2nd or- 
der analysis. Only a monotonic enhancement of the scaling 
properties is observed. 

In Table |3| we compare the results from different approx- 
imations to the exact dynamics in the perturbative regime 
(exact PT) , available in the literature and we give as well the 
SC predictions in the same regime. We recall that the FFA 



Dynamics 


S2,4 


•S'3,0 


53,2 


•S'4,0 


54,2 


FFA* 


0.43 


3 


1 


16 


15.0 


LPA* 


0.72 


3.40 


2.12 


21.22 


37.12 


ZA* 


1.27 


4 


4.69 


30.22 


98.51 


SC 


1.44 


4.86 


10.08 


45.89 


267.72 


Exact PT* 


1.82 


4.86 


9.80 


45.89 


7 



Table 3. Comparison between different non-linear approxima- 
tions to the for the unsmoothed fields up to the first corrective 
term beyond tree-level (the one-loop term). The asterisk denotes 
the results obtained within the diagrammatic approach for the 
relevant dynamics. 



is based on a linearization of the peculiar velocity field. The 
LPA assumes the potential remains linear throughout the 
gravitational evolution, and the ZA takes the trajectories of 
the particles to be straight lines. All values concerning the 
one- loop contribution (including the exact PT ones), except 
for the SC ones, are derived making use of the diagrammatic 
formalism and are given in SF96a. Tree-level amplitudes had 
been derived previously though and a summary of those re- 
sults may be found in Bernardeau et al. (1994). As summa- 
rized in Table 3, the SC model yields the best estimates for 
the loop corrections in PT within the available non-linear 
approximations. 



5.2 Previrialization in the Variance 

It is important to notice the different behavior between the 
cumulants (like the variance) and the hierarchical ampli- 
tudes, Sj. It has long been argued that non-local evolution- 
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10 

R (Mpc/h) 



50 



Figure 3. The one-loop contribution to the variance in PT. 
The sohd hne depicts the exact PT behavior, the short-dashed 
hne shows the shearless (SC) contribution while the long-dashed 
line gives the tidal contribution. Two characteritic indices arise 
in this analysis: ric — —1.4 where the one-loop variance vanishes, 
and n* ~ —1.6, where tidal effects are zero and the SC model 
prediction is exact. 



Figure 4. The variance in spherical shells of radius R. We com- 
pare linear theory to the SC model and the exact PT (leading 
order) non-linear correction, against the linear term. Each pan- 
els, (a) and (b), show the ratio of the non-linear to the linear 
variance for the P = 0.2 (a) or P = 0.5 (b) CDM initial P(k). In 
both cases the continuous line from bottom to top corresponds 
to (Tg = 0.4, 0.6, 1.0. The SC predictions for the same amplitudes 
are shown as dashed lines. 



ary effects generate a suppression of collapse on the largest 
scales due to the increase in small-scale power that generates 
large scale random motions, the so-called effect of previri- 
alization (see Davis and Peebles 1977, Peebles 1990, Lokas 
et al. 1996). This will result in a non-linear variance that is 
smaller than the linear one within a certain range. This ef- 
fect has been found in N-body simulations {e.g., see Figure 
9 inBGE95). 

The local picture arising from the SC model is unable to 
account for a complete suppression of non-linearities in the 
large scale density (and velocity) fields and yields a faster 
growth of fluctuations on the large scales compared to the 
actual non-local dynamics. This qualitative feature of the SC 
model was already pointed out by Bernardeau (1994c), in 
what he called the rare event limit, 5i/ai -^ cx). Making use 
of the third order in PT (in Lagrangian space) for the density 
field smoothed with a top-hat window (see his Eq. [38]) for a 
scale free power spectrum, he claimed that the SC picture 
deviates from the exact PT as smoothing effects increase. In 
particular, he found that for a spectral index n > —1, PT 
diverges due to the shear effects. Whenever more realistic 
spectra (models with a cut-off on small scales, such as CDM 
or APM-like ones) are introduced, this divergence becomes 
only logarithmic. 

In PT, the above commented suppression depends on 
the smoothing effects and leads to the appearance of a crit- 
ical value of the spectral index (assuming scale free power 
spectrum) around Uc ~ —1.4 (see S97), for which non- linear 
contributions vanish. Figs |l| and y illustrate well this point 
since no change of sign in the first non-linear contribution 
(beyond tree-level) to the evolved variance is appreciated in 
the plot contrary to what has been found in numerical sim- 



ulations depicting the exact perturbative scaling (see Lokas 
et al. 1996, Fig 4). Note that the SC prediction yields a 
minimum contribution at n ~ —1.14, which is not far from 
ric ~ —1.4, but the net effect at the minimum is non-zero. 

Comparing the results for the variance from the ex- 
act PT and the shearless (SC) approximation, we can in- 
fer some information about where the tidal effects become 
lower, comparable or greater than the shearless contribution, 
provided the variance changes smoothly. Let us formally de- 
compose the one-loop contribution to the variance in the 
shearless (SC) and tidal contributions as follows. 



^2,4 



SC , Tide 
^2,4 + S2 4 , 



(57) 



where sf 4 > for all n as displayed on the bottom panel of 
Fig hi The one-loop contribution to the variance in PT, s^J 
can be integrated from the analytic expression for the one- 
loop power-law power spectrum given in SF96b (see their 
Eq.[B5]) valid in the range — 1 > n > —3. As shown in Fig. 
bl the shearless (SC) contribution gives a systematic un- 
derestimation of the PT value from the unsmoothed value 
(n — —3) up to n* ~ —1.6 where the tidal contribution 
vanishes, and thus, s^J — sf 4 . Around this index, the SC 
model gives an accurate estimation of the one-loop variance. 
Beyond this point, n > Ui,, the tidal contribution rapidly 
grows and becomes negative (see Fig.0). For Uc — —1.4, 
sf 4 — —82^4" and the one- loop variance vanishes. As seen 
in the plot, while the shearless contribution remains positive 
and small, the tidal contribution takes larger and larger neg- 
ative values as n ^ —1, what breaks down the PT approach 
(at least for the variance). 

In Figure H we compare the non-linear evolution of the 
variance for the SC to that of the exact PT. We present 
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the ratio of the non-linear to the linear variance from CDM 
model, with T = 0.2 (top panel) and T = 0.5 (bottom), 
as a function of the smoothing radius. The leading or- 
der non-linear SC predictions are shown as dashed lines. 
The corresponding (numerically integrated) PT results are 
shown as continuous lines. These are estimated from P{k) 
to second order (2nd non-vanishing contribution) as ob- 
tained following Baugh & Efstathiou (1994). In both cases 
we show results for different amplitudes of the linear vari- 
ance: (78 = 0.4,0.6, 1.0, which has a linear variance a ~ 0.5 
at scales 7? ~ 6, 10, 15 h~^ Mpc, respectively. We can define 
an effective index, n^ff, as the index at the scale where the 
variance becomes unity. For the F = 0.2 model we have that 
n^ff'^ — —2, while for F = 0.5 this is about 0.5 larger (this 
can be seen for example in Figure 9 in Croft & Gaztaiiaga 
1998): nlff-^ ~ -1.5. The SC prediction should match bet- 
ter the F = 0.2 model, as the effective index is closer to 
n*, where tidal effects vanish. For the F — 0.5 model we 
have that Usff — ric and the SC model is not such a good 
approximation. In this case we expect non-linear effects to 
be small due to the cancellation of tidal and purely local 
contributions to the cumulants. This can be seen in Figure 
kl The SC model matches well the PT prediction for the 
F = 0.2 model, where non-linear effects are more impor- 
tant. It also matches well the earlier evolution of F = 0.5 
model, for small amplitudes of oi, e.g., earlier outputs given 
by smaller a^. For the later outputs {e.g., as ~ 1 in Figure 
kl bottom panel) tidal effects become important canceling 
out the non- linear growth (previrialization effect). This is 
not fully recovered by the SC model (which only accounts 
for the local effects). Nevertheless note that for small scales, 
where the variance is large, the prediction seems to be dom- 
inated by the local effects and the SC model also recovers 
the exact prediction. 

In Figure ^ we show the non-linear evolution of the 
variance as the ratio to the linear variance from CDM N- 
body simulations with F = 0.5 (filled triangles) and F = 0.2 
(open triangles) as a function of the smoothing radius. The 
error-bars, from 10 realizations of the models, which are not 
displayed for clarity, are always smaller than the symbols 
except for the last 3 points where the error-bars are smaller 
than twice the size of the symbols. In both cases erg = 0.5 
and the linear variance CT~lati?~2/i~^ Mpc. The SC 
predictions are shown as dashed lines, which are to be com- 
pared to the exact PT results (continuous lines). For the 
F = 0.2 model there is hardly any difference between the 
predictions and the N-body results. In fact, the latter fol- 
lows both predictions up to large values of a (small scales) . 
For the F = 0.5 model, non-linear effects are small during 
the weakly non-linear phase and both linear and non-linear 
predictions are close to the N-body results. There is a slight 
sign of previrialization at _R ~ 20/i~^Mpc, but the effect 
is quite small. This way, despite the exact PT results are 
more accurate in qualitative terms, there is little difference 
in practice. Note that, in this case, the agreement with pre- 
dictions does not extend to large values of a (small scales) 
unlike the case of the F = 0.2 model. 

Figure shows the corresponding results for o-g = 1. 
In this case, the F = 0.2 CDM model is the one with an 
effective index rie// — Wc so that non-linear effects are small 
and we have a similar situation to the F = 0.5 model at 
erg = 0.5. The effective index for F = 0.5 model at erg = 1.0 




R (Mpc/h) 

Figure 5. Non-linear evolution of the variance. Symbols show the 
ratio of the non-linear to the linear variance from CDM N-body 
simulations with F = 0.2 (top panel) and F = 0.5 (bottom panel) 
as a function of the smoothing radius. In both cases erg = 0.5. The 
SC model predictions are shown as a short-dashed line while PT 
predictions are shown as a continuous line. The arrows indicate 
where ct; = 0.5. 




Figure 6. Same as Fig. H for erg = 1.0. 



si n > nc so that the previrialization effect discussed above 
is larger and the SC prediction fails (long-dashed line). The 
prediction for the second order contribution to P(fc) for 
F = 0.5 is shown as a continuous line. Although it shows 
some previrialization effect (the ratio is smaller than unity) , 
it only fits the N-body results in the narrow range of scales, 
R > 30 /i~^ Mpc, where linear theory is a good approxima- 
tion, given the errors. The above arguments explain in a sim- 
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Figure 7. The variance for APMPK at crs = 0.8 (figures) com- 
pared to PT models (lines): (a) linear theory prediction (continu- 
ous line), (b) the 2nd order result, i.e., when the 2nd perturbative 
contribution to the variance is included in the SC (short dashed 
fines) (c) the 3d order result in the SC (long dashed lines). 



pie fashion why the SC model matches so well the variance of 
the APM simulations even on small scales (for output times 



0-8 



< 



1) as shown in Fig. M. The effective index rie// on the 



> 



-1.5, a 



quasi-linear scales is within the range — 2 > n^ff 
bit lower than the T = 0.2 CDM model (see Fig.3 in BG96, 
where quasi-linear scales correspond to k ~ 0.1). For this 
index, tidal effects are expected to give a sub-dominant con- 
tribution. Although the SC model seems to match well the 
simulations up to the scales where the linear variance cr; > 1, 
e.g., beyond the point where PT must break down, the third 
order SC correction (long-dashed line) shows significant de- 
viations already at J? ~ 6 h~^ Mpc, where a; ~ 1. 

The results here obtained are in qualitative agreement 
with those found by B94c, according to which, the rare event 
limit progressively differs from the exact PT predictions as 
n — > —1, even in the first correction to the variance S2,4. 
Furthermore, B94c also finds no indication of previrialization 
in the rare event limit and stresses the unimportance of this 
effect for the APM & IRAS observations, given the effective 
spectral index related to the scales on which they lie {rieff — 
— 1.3,-1.4, respectively). However, according to B94c, the 
rare event limit is always different from the exact PT result 
whenever Si ^ 0. This is at variance with what we find here, 
where the SC model is expected to yield accurate results, 
particularly for values of the spectral index n £ [—1.5, —2]. 
What is more, there must be some value n* within that 
interval for which tidal effects exactly vanish. 



5.3 Comparison of the Sj Ratios with N-body 
Simulations: Previriahzation Lost 

On the other hand, when it comes to evaluating the hier- 
archical amplitudes, Sj, it turns out that the cancellation 
of (non-local) tidal effects erases the previrialization effect 
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Figure 8. The hierarchical skewness 53 from 10 realizations of 
flat CDM N-body simulations at output time eg = 0.5 (top panel) 
and erg = 0.7 (bottom panel). Symbols with error-bars correspond 
to r = 0.2 (squares), and P = 0.5 (triangles). Each case is com- 
pared to the corresponding PT predictions (lines): a) the PT tree- 
level prediction uncorrected for the ZA transients (dot-dashed 
line), and corrected for the transients (continuous), b) the 2nd 
order result in the SC uncorrected for the ZA transients (short- 
dashed), corrected for the ZA transients (long-dashed), and the 
case with 72 = (dotted). The arrows show where a = 0.5. The 
dot-dashed line shows the tree-level prediction for the P = 0.5 
and erg = 0.7, uncorrected for the ZA transient. 



and the suppression of non-linearities at the first corrective 
order beyond tree-level appears at a different spectral in- 
dex n « 0, what seems to be an intrinsic (local) shearless 
feature if the good agreement between the local (SC) and 
non-local (exact PT) hierarchical amplitudes is anything to 
go by. Furthermore, loop-calculations in Fourier space con- 
cerning the reduced Bispectrum, Q, seem to confirm this ar- 
gument as the scale dependence of Q retains non-vanishing 
non-linearities even at n ~ —1.5 (at least for the equilateral 
configuration, see Scoccimarro et al. 1998, Eq.[30]). How- 
ever, the reduced bispectrum also preserves some non-local 
features which manifest in terms of a loss of configuration 
dependence (isotropization) roughly at the same value of the 
spectral index as the one estimated from the variance or the 
power spectrum, Uc ~ —1.4 (see SF96b, and S97). 

In Figs y, we compare the SC predictions with CDM 
simulations for two different models. The smoothing coeffi- 
cients Ffe are modified to include the h ighe r- order logarith- 
mic derivatives of the variance (see Eq. ]A3| in Appendix Mj 
but it turns out that, within the errors, they do not signif- 
icantly change the predictions derived from the power-law 
power spectrum. At small scales 72 changes 5*3 by less than 
10%, which hardly makes any difference given the errors, as 
shown by the dotted line, only displayed for the F = 0.2 
erg = 0.05 model for clarity. The error bars shown in the 
SCDM plot are derived from 10 realizations. The arrows in 
the Figure show the scale at which a = 0.5. 

On large scales the spectral index n for CDM models 



© 0000 RAS, MNRAS 000, 000-000 



16 P.Fosalha and E.Gaztanaga 



becomes larger and non-linearities are very small, as pre- 
dicted by the SC model (see Table m). This explains why 
deviations from the tree-level are hardly noticeable in the 
later outputs (bottom panel) and they only become impor- 
tant, given the errors, when PT theory should break down 
[ai ~ 1). The effective index in the earlier outputs is smaller 
(corresponding to smaller scales) and deviations from tree- 
level are more significant. In this regime the SC predictions 
show a very good agreement with the N-body simulations 
(top panel) up to (t ~ f. Note that the predictions for the 
exact PT corrections are not available in this case, but the 
agreement with the N-body simulations indicates that the 
SC models provides an excellent approximation. 

We have corrected for ZA transients, an artifact due to 
the initial (ZA) start in the N-body simulations (see BGE95 
1995, Scoccimarro 1998). We use the analytic results of Scoc- 
cimarro (1998) for the tree-level: 



5^3,0 = y + 7 - 



16 
35' 



-7/2 



(58) 



where a is the expansion factor away from the (ZA) IC. In 
principle this is not a very important correction for the tree- 
level result, as the sampling errors on large scales, where a 
is small, are large. This is more important for smaller scales, 
as the errors are smaller. In Figs H we show the tree- level 
prediction for a ^ oo (dot-dashed line) as compared to the 
actual prediction, including the transient, for a — 5.0 (the 
bottom continuous line) in the F — 0.5 model for erg = 0.7. 
The long-dashed line shows the 2nd order (one-loop) pre- 
diction in the SC model corrected for the ZA transients. 
This can be done by using the tree-level results for the tran- 
sients, as they determine all higher-orders in the SC model. 
As pointed out by Scoccimarro, the effects of the transients 
tend to create the false impression that the tree-level PT 
has a wider range of validity. In our case, by comparing the 
points to the dot-dashed line one might conclude that tree- 
level results are valid up to i? ~ 9 h~^ Mpc, while a compar- 
ison to the corrected prediction (continuous line) indicates 
i? ~ 15 h~^ Mpc. On the other hand, adding the loop terms 
to the uncorrected tree-level would give the false impression 
of a poor agreement with the SC model. 

Figs y and hd display the scaling of the higher-order 
moments according to the APM-like simulations (the mean 
of APMPKl and APMPK2 in Table 0). As for the previous 
cases, the SC model matches well the observed behavior for 
the higher-orders in the N-body simulations up to the scales 
where PT is expected to break down, e.g., ai « 1. As we 
consider higher order moments, deviations from the tree- 
level prediction are larger, a trend that is reproduced by the 
SC predictions. We have found a similar trend for higher 
orders in the CDM models. 

The ZA transients, shown as dot-dashed lines in Fig. bl 
are not very important in this case (a ~ 6). The effect of 
using 72 = is also small although it tend to introduce a 
slight shift in the predictions. 

These results extend the domain of applicability of the 
SC model in PT beyond the point expected from B94c anal- 
ysis of the rare event limit, as potential cancellations of the 
tidal contributions in the Sj ratios were not examined there. 




6 8 10 20 

R (Mpc h"') 



Figure 9. The skewness, 53, and kurtosis, 54 for simulations of 
the APMPK at ag = 0.8 (open squares with error-bars) compared 
to tree-level theory prediction (continuous line) and the 2nd order 
result {Sj including the one-loop term) in the SC (short dashed 
lines). The corresponding predictions for 72 = are shown as 
long-dashed and dotted lines respectively. The dot-dashed lines 
show the tree-level predictions uncorrected for the ZA transients. 



6 DISCUSSION AND CONCLUSIONS 

No available analytic approximations to the dynamics of cos- 
mological perturbations lead to accurate predictions for the 
one-point cumulants in the non-linear regime (see §3.7). This 
applies to all popular approximations, such as those based 
on putting some constraints on any of the fields that couple 
to the density (such as the FFA or the LPA), the Zel'dovich 
Approximation (ZA), which puts a constraint on the trajec- 
tories of the masses, or the local Lagrangian approximations, 
which assume that the evolution of the density contrast may 
be described locally, i.e., without any infiuence of the sur- 
rounding matter to the mass trajectory. Despite giving some 
useful insight to the exact dynamical picture with much sim- 
pler calculations, none of them is able to reproduce the exact 
values of the one-point cumulants even at tree-level as de- 
rived in the exact PT. 

As mentioned in the introduction, there are serious 
difficulties in applying the exact PT approach to find the 
gravitational evolution of the cumulants of a Gaussian field 
to higher-order (loop) corrections, i.e., to compute next- 
to-leading orders. For non-Gaussian initial conditions this 
problem is apparent even to leading order (in non-linear cor- 
rections). An important simplification is found when only 
the monopole or spherical contribution is considered. This 
contribution is exact for all the tree-graphs. We have shown 
how the solution in this case can be found directly from 
the dynamical equations of the evol ution of a cosmic fluc- 
tuating field (see e.g., Appendix 4.3 and Appendix O) as it 
is given by the spherical collapse (SC) model for the case 
of density fiuctuations for a non-relativistic presureless ir- 
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Figure 10. Same as Fig. H for Ss and Se. 



rotational fluid. This provides us with a simpler derivation 
and interpretation of the results presented by Bernardeau 
(1992), showing why the SC model gives the exact tree- level 
in PT. This lacked a satisfactory explanation in the work of 
Bernardeau. 

We have explored the predictions for the one-point cu- 
mulants of the density field in the weakly-nonlinear regime 
within the spherically symmetric approximation to the dy- 
namics. This is done both for the exact dynamics and for 
the Zel'dovich approximation (see Appendix O). As the SC 
model gives the exact tree-level irrespective of the nature of 
the initial fluctuations and the value of Q, it can also be used 
to make predictions for non-Gaussian initial conditions (the 
results are presented in the accompanying Paper II), the ve- 
locity field and the f2-dependence of these predictions (see 
Paper III of the series). 

A natural shortcoming of the SC model is the loss of 
previrialization, i.e., the generation of large scale random 
motions induced by the small scale power, which is a non- 
local phenomenon. This is in qualitative agreement with the 
analysis given in B94c in the rare event limit to PT which 
is described by the SC model. The lack of previrialization is 
illustrated by the absence of vanishing non-linearities in the 
variance for a critical index Uc ~ —1.4 (see §[3.2|), which is 
a direct consequence of previrialization. This phenomenon 
is expected from the exact PT (SF96b) and it has been ob- 
served in numerical simulations (Lokas et al. 1995). 

Despite this, we have shown that the SC model also 
yields accurate predictions for the cumulants beyond the 
tree-level. Tidal effects in the variance are actually small 
because the effective index n^//, as measured in galaxy cat- 
alogues (such as the APM), is n^ff — [—1,-2] (see §5.2) 
for which either non-linearities are small or tidal effects are 
sub- dominant. 

We stress the importance of applying the SC approx- 
imation in Lagrangian space. There, the SC model is de- 



scribed by a transformation that only depends on the value 
of the linear field at the same point (what we call a local- 
density transformation). However, when going back to Euler 
space the density fluctuation (defined at a point) is normal- 
ized with a factor (see Eq.|p9|) which is a function of the 
(non-linear) variance. Since the variance is a volume average 
of the two-point correlation function, this factor yields some 
non-local contribution to the cumulants (in Euler space). 
This non-local contribution is missed when introducing the 
SC model in Euler space directly, thus is not surprising that 
the predictions for the cumulants in the SC approximation 
in Euler space are a poor estimation of those in exact PT, as 
the latter are dominated by the non-local (tidal) effects (see 
Table H in Appendix pl). The contribution to the cumulants 
of this non-local term is typically negative so the predictions 
from the SC model in Euler space generically overestimate 
those in Lagrangian space. 

For the predictions within the SC model (in Lagrangian 
space) for the hierarchical ratios, Sj, tidal effects partially 
cancel out, at least to one-loop. We find an excellent agree- 
ment for the hierarchical amplitudes Sj in the perturbative 
regime of the SC model, with those derived by SF96a for 
the exact PT in the diagrammatic approach. A similar con- 
clusion follows for the Zel'dovich dynamics (Appendix pi), 
where the monopole approximation (the SSZA) results are 
in excellent agreement with the exact calculations. 

We have also compared the predictions for the higher- 
order moments from the SC model to those measured in 
CDM and APM-like N-body simulations, and they turned 
out to be in very good agreement in all cases up to the scales 
where ai ~ 1, supporting our view that the tidal effects 
have only a marginal contribution to the reduced cumulants. 
Furthermore, the break down of the shearless approximation 
roughly coincides with the regime for which the perturbative 
approach itself breaks down, a ~ 0.5. That is, where the 
contribution of the second (one loop) and the third (two 
loop) perturbative order in the SC model are significantly 
different. 

These results extend the domain of applicability of the 
SC model in PT beyond the point expected from B94c anal- 
ysis of the rare event limit, as potential cancellations of the 
tidal contributions in the Sj ratios were not examined there. 

Future work intended to extend the SC model should 
try to incorporate the local contribution to the tidal field by 
including other multipoles in the kernels _F„, e.g., including 
the shear by using the dipole contribution. 
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APPENDIX A: CUMULANTS FROM THE 
LOCAL-DENSITY TRANSFORMATION OF 
GAUSSIAN INITIAL CONDITIONS IN 
LAGRANGIAN SPACE 

For a generic model of non-linear evolution of the density 
field in which all the information is encoded in the linear 
density field alone, one can construct a transformation of 
the kind. 



^[^^] - E ?} [^'] 



(Al) 



Scoccimarro, R., Frieman, J., 1996a, ApJS, 105, 37 (SF96a) 
Scoccimarro, R., Frieman, J., 1996b, ApJ, 473, 620 (SF96b) 
Scoccimarro, R., Colombi, S., Fry, J.N., Frieman, J., Hivon, 

E., Melott, A., 1998, 496, 586 
Zel'dovich, Ya.B., 1970. A & A, 5, 84 



(see Eq.fla]) with ci = 1, to reproduce the linear solution. 
This is what we shall call the local- density transformation. If 
we further assume that the above transformation describes 
the non-linear evolution in Lagrangian space, it has to be 
normalized 

Next we present the general results for the cumulants 
including their first non-vanishing perturbative orders for 
GIC. For the SC model the calculation corresponds to the 
process described in M Note that this can be extended in 
a straightforward way to NGIC, just by taking into account 
the relevant terms to be kept in the perturbation series (see 
Paper II). 

We proceed and give the coefficients of the perturbative 
terms according to the notation introduced in 95l 

S2,4 = 3 - 4C2 + C2/2 + C3 

S2,6 = 15 - 36c2 -I- 81ci/4 - 3ci/2 + lOca - 7c2C3 

+ 5c^/12 - 7C4/4 + C2C4/2 + C5/4 

S-ifi = 3c2 

S3,2 = (-4+15ci-4ci-8c3-|-3c4)/2 

53,4 = -18 + 44C2 -I6C2 -3C2 -7C2+5C2/4 

- 22C3 + 12C2C3 + 10C^C3 + S,C-A - Acj 

- 7C2C3/4 + 5C4+9C2C4/8-9C2C4/4 
+ C3C4 — 2C5 -|-3c2C5/4 + 3C6/8 

^4,0 = I2cl + 4C3 

S'4,2 = 6 - 36c2 4- 15ci + 60ci - 15c2 - 4c3 

- 2OC2C3 — 6C2C3 — 5C4 + I8C2C4 + 2C5 

54,4 = (180 - 984c2 + 1488ci - 534ci + 99c2 

- I8OC2 -I- 27c2 + 260c3 - 8I6C2C3 

+ 232cic3 + 144cic3-H30c|c3-H76c^ 

- 48C2C3 - 30cic3 - 4C3 - 86C4 

-|- 163c2C4 + I29C2C4 — 66C2C4 — 72C3C4 

-f 8C2C3C4 -I- I6C4 + 8C5 — 4OC2C5 
+ I2C2C5+7C3C5 -5C6 + 12C2C6+C7)/2 
55,0 — 6OC2 + 6OC2C3 -I- 5C4 

85,2 = -24 + 180c2 - 5IOC2 4- 240c2 -I- 450c2 

- 108c^ - 8OC3 + 9OC2C3 - 120cic3 

- 9OC3 - 10c4 - 105C2C4/2 + 170cic4 
+ 5OC3C4 - 6c5 + 4OC2C5 -I- 5c6/2 

Ss.o = 6(604 + 120cic3-hl5c^-h20c2C4 + C5) 

Sfi,2 = 120 - IO8OC2 + 3870C2 - 6930C2 + 3060c2 

+ 3600C2 -840c2 +480C3 -3240C2C3 

+ I5OOC2C3 + 3840C2C3 - I8OOC2C3 
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IIOC3 - I62OC2C3 - 225cic3 + 304 
I5OC4 - 135C2C4 - 6OC2C4 + 1500cic4 
42OC3C4 + I5OOC2C3C4 + 105c| - I8C5 
IO8C2C5 + 585c2C5 + I35C3C5 — 7c(i 

75C2C6 + 3C7. 



(A2) 



Note that all the values of Cn can be obtained (or rewritten) 
in terms of the tree-level alone, Sj,o, indicating that the 
knowledge of the tree-level is enough to fully specify the 
underlying local-density transformation, and, therefore, to 
generate all higher order corrections. 

Furthermore, we can make use of the general formula 
for the top-hat filtering derived in §4.4, and generalize the 



above given expressions by replacing the unsmoothed Ck co- 
efficients by their smoothed counterparts Cfc in the following 
way: 



C2 



C3 



C4 



C2 



Ci 



21 

3 

y + -C2 7l + - + y 



-I- 



4 . „ 2 8 

C4 + - 71 - 4 C2 71 + 2 C2 71 + - C3 71 

4 2_^8 2_^ 8 3 -^ 

3 71 + 3^271 + ^71 

4 2 

g 7172 + ^73, 



3 72 + 3 C2 72 



(A3) 



and so on, where: 



7p = 7p(-R) 



d logf' R ' 



These expressions are valid for arbitrary dynamics (to be 
specified through the unsmoothed coefficients) and a generic 
initial power spectrum. They can be used for the particular 
cases of interest. For the SC dynamics {ck = Vk as given by 
Eq.M]) we find for the first coefficients of the cumulants. 



S2,4 = 
^3,0 = 
53,2 = 

+ 
54,0 ~ 
5*4,2 = 

+ 
+ 



1909 143 11 2 , 72 

7l H 7i + -^ 

1323 126 '36 '^ 6 
34 

y + 71 

1026488 12862 407 2 10 3 

101871 "^ 1323 '^^ '^ 126 '^^ "^ 27 "^^ 

11 _L 2 73 

y 72+ 37172 + y 

60712 62 7 2 2 

3^ + ^71 + 371 + 372 

22336534498 42649448 3571621 2 

+ -^77:^::^ 7i + 0001. Ti 



83432349 130977 

35047 3 1549 4 . 575777 



1134 
263 2 



7i + 



648 



7i + 



25 



5^7172 + ^72 + 



11907 
2084 86 



23814 

5981 

72 + 



189 



7172 



— 73 + ^7173 + ^74 



(A4) 



This way, the results for the SC and SSZA dynamics de- 
scribed in the text (see §0 and Appendix O respectively) are 
given as particular cases of this local- density transformation 
in Lagrangian space whenever we replace the cj, coefficients 
by those associated to the relevant dynamics: Ck — Uk or 
Ck =T^ for the unsmoothed or smoothed fields respectively 
(the previous being a particular case of the latter) . The same 
expressions hold for the velocity divergence fields, replacing 



Ck = Uk OT Ck = jMk for the unsmoothed and smoothed fields 
respectively. These fik coefficients are related to those for the 
density field through the equation of continuity (see Paper 
III for details). 



APPENDIX B: CUMULANTS FROM THE 
LOCAL-DENSITY TRANSFORMATION OF 
GAUSSIAN INITIAL CONDITIONS IN EULER 
SPACE 

In this section we derive the cumulants from the local- 
density transformation of GIC in Euler space to quantify 
the departures form the Lagrangian formulation and give 
values for the SC approximation. The formulae given below 
reproduce those given in Fry & Gaztafiaga (1993). There, 
they were presented as a bias transformation between the 
luminous and the underlying matter fluctuations. They can 
also be obtained by simply replacing _F„ = Cn/n\ in the ex- 
pressions for the loop corrections (in Euler space) given in 
SF96a. 

The first perturbative contributions to the cumulants 
in Euler space are the following. 



S2,4 = 


- | + « 


53,6 = 




= cl+ 6C2C3 + -C4 


53,0 = 


= 3C2 


53,2 = 


= -2cl + |c4 


54,0 


= 12C2 + 4C3 


54,2 = 


— 15C2 — 6C2C3 + 



IC4 + 2C5. (Bl) 

where we have also computed the third order cumulant, ^3 
whose one-loop contribution is denoted by ^3,6 following the 
notation introduced in ipl It is easy to see that this contri- 
bution is given by ^3,6 = 53,2 + 2s2,453,o- 

Replacing in these expressions the values for the SC dy- 
namics with a top-hat filter and a power- law spectrum, i.e., 
introducing the smoothed coefficients given by Eqs.[k7|,[B7|, 
we get estimates for the cumulants as summarized in Table 



Bl. For a comparison of both the Lagrangian and Eulerian 
predictions within the SC model with respect to the exact 
FT values we also display the results by SF96a, SF96b and 
S97 obtained using the diagrammatic representation for PT. 
Notice that unlike the case for the SC estimates from 
Lagrangian space, the Euler estimates just give the right or- 
der of magnitude for the cumulants when compared against 
exact analytic calculations (see Eqs.[|53[, |55|). The observed 
trend is that (non-reduced) cumulants in Euler space typi- 
cally overestimate those in Lagrangian space as the normal- 
ization factor (see Eq.M]) between the two spaces intro- 
duces terms that give a negative non-local net contribution 
in terms of the variance (see §5.1 for a discussion). 



Since this overestimation effect in the higher-order cu- 
mulants scales less steeply than hierarchically with respect 
to that of the variance, A^j < {AS,2)''~^, the effect in the 
Sj ratios is thus dominated by that in the variance. This 
results in a net underestimate of the Sj ratios. Notice that 
the relative underestimation becomes larger and larger as 
smoothing effects increase. 
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Smoothed 






7 = 


7 = -l 


7 = -2 


7 = -3 




n = —3 


n = -2 


n = -1 


71 = 


■^2,4 


4.92 


2.76 


1.20 


0.26 


■^2,4 


1.44 


0.61 


0.40 


0.79 


*2,4 


1.82 


0.88 


7 


7 


/zE 


54.64 


21.80 


5.68 


0.38 




24.09 


7.95 


2.85 


2.91 


cPT 
?3,6 


27.49 


9.97 


7 


7 


"^3,2 


6.85 


0.54 


-1.21 


-0.60 


•^3,2 


10.08 


3.21 


0.59 


-0.02 


qPT 

"^3,2 


9.80 


3.18 


7 


7 


cB 

•^4,2 


208.45 


24.83 


-10.71 


-3.09 


•^4,2 


267.72 


63.56 


7.39 


-0.16 


cPT 
•^4,2 


7 


? 


7 


? 



Table Bl. Comparison between the cumulants in the SC model 
in Euler space (^) with those in Lagrangian space (^) and in 
exact PT {^^). 



APPENDIX C: TOP-HAT SMOOTHING IN 
FOURIER SPACE 

To illustrate the fact that the SC model gives the correct 
smoothed tree-level amplitudes for Gaussian initial condi- 
tions for a top-hat window, we shall derive TTJ explicitly 
by imposing the spherically symmetric approximation in 
Fourier space and show that the skewness S3 = 37^, ex- 
actlyreproduces the leading order exact perturbative result, 
Eq.M]; e.g., P2 = i/2 + 7/3 We perform the calculation in 
Euler space but the result does not differ from that in La- 
grangian space at tree-level for Gaussian initial conditions, 
as discussed in 94. Differences are only expected to appear 
in higher perturbative orders (cr-corrections). 

To see this, we recall the properties derived by 
Bernardeau 1994b (B94b henceforth) for a top-hat window 
function (spherical window) defined as, 

Wth{X) = 1 if |x|<i?o, 

and otherwise, for a scale Ro, so that the smoothed fields 
are obtained through the integral, 

5 = S{Rq) = / cPx Wth{x)5{X), 

We turn to Fourier space for convenience, where the 
smoothed fields are expressed as. 



5{Ro) = 



d'k 



WTH{kRo)S{k), 



(27r)3/2 
where WrHikRo) is the Fourier transform of Wth{X). In 



particular, for the second-order in the perturbative series, 
we have for the smoothed field 



52iRo) = 

where, 
Pi,2 = 1 



d'^ki cPk2 



(27r)3/2 (27r)3/2 

fcl ■ K>2 



5k,5k,_W{\k^ + k2\Ro) 



\Q1.2) + \d2Qi,2 



kl 



= 1 + 



(fcl ■ k2f 



and Di are the time-dependent growth factors for the i-th 
perturbative order to be solved with the SC equations of 
motion. Now, decomposing the integrals into their angular 
and radial part and translating the property of spherical 
symmetry into Fourier space language, 



5k. = S, 



fed' 



(CI) 



we can apply some properties for the top-hat window func- 
tion (see B94b, Eqs.[A5] and [A6]) and get. 



52{Ro) 



5i{Ro) = 



O2 1 d \og[5i{Ro)Y 

2Dl 6 dlogi^o 



MRo)f 



d'^k 

(27r)3/: 



WTH{kRo)5i:Di{t). 



We further use the fact that in general the smoothed linear 
density contrast at a point x, can be factorized in its scale 
dependent part a{Ro) and its normalized linear field e{x), 
5i{Ro) = 5i{Ro,x) = ai{Ro) 'ri{^)j ^nd we finally get: 



52iRo) 



34 1 d\og[ai{Ro) 
21 3 d log Ro 



[Si{Ro)f 



If we define, (52 (i?o) 



-u^/2) [Si{Ro)Y , we and 

2 



_ 34 1 dlog[cr,(it:o)]- , 7 

21 3 dlog_Ro 3 



(C2) 



(C3) 



which exactly reproduces Eq.pTf. This result can be ex- 
tended to higher orders following the properties of the top- 
hat window presented in B94b. 



APPENDIX D: THE SPHERICALLY 
SYMMETRIC ZEL'DOVICH APPROXIMATION 

A simple accurate description of the non-linear dynamics 
of the fluid elements before shell crossing (single streaming 
regime) is provided by the so-called Zel'dovich approxima- 
tion (hereafter ZA, see Zel'dovich 1970). According to this 
approximation, particle positions in comoving coordinates 
are assumed to follow straight lines. Despite being exact 
only at linear order in Lagrangian FT, the latter solution 
has been successfully tested as an accurate approximation 
for the description of the dynamics in the non-linear regime 
(see e.g., Croft & Gaztanaga 1998). Within the ZA the fluid 
equations can then be easily integrated, and yield 



i+.=n 



(l-D(i)AO' 



(Dl) 
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where N is the number of spatial dimensions and Xi = 
—dipi/dqi, which is in general a non-local relation between 
the displacement field and the evolved density fluctuation. 
At this point we introduce the spherically symmetric as- 
sumption for the ZA dynamics (SSZA) as we did before with 
the exact dynamics. Thus we set Ai = A, i.e., all directional 
derivatives are equal due to the spherical symmetry (no tidal 
forces) and rewrite, 5i = D{t) ^ ,_ Ai = 3 D{t) A, as follows 
from linearising Eq. [Dl]. The assu mption of spherical sym- 
metry renders the relation Eq. |D1[ as a local transformation 
that may be expanded in power series to get the full pertur- 
bative series for the density contrast. In particular, for the 
3D case, the normalized SSZA transformation has the form. 



Sr^ii-Si/sy 



(D2) 



to be normalized according to Eq.|39|]. Notice that the latter 
transformation may be straightforwardly extended to any 
spatial dimensions. As seen in Fosalba et al. (1998), the ex- 
ponential transformation that leads to the lognormal hier- 
archical amplitudes, Sj = J^~^, are only recovered when 
one takes the limit of infinite spatial dimensions (see also 
Bernardeau and Kofman 1995), which does not make much 
physical sense. 



Expanding Eq.|D2 in power series we can determine 
the unsmoothed coefficients v^ of the local non-linear trans- 
formation for the SSZA, 



V2 = -^ 1.33; 
^•4 = — ~ 4.44; 



;/3 = ^ ~ 2.22 
9 



I'b 



280 
27 



10.37, 



(D3) 



and so forth. Introducing them in the smoothed transforma- 
tion according toEq.pTf for a top-hat window, we can now 
find, using Eq.|A3 with c^ = Vk, the following results for 
the smoothed density field for GIC, and a power-law power 
spectrum. 
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The results for the variance and S,j for the SSZA fo r pa rtii 



Dl 



ular values of the spectral index are given in Table 

For the unsmoothed fields (7 = 0), analytic results in- 
cluding the two corrective terms beyond the tree-level, were 
derived by SF96a in the context of the diagrammatic ap- 
proach. For the variance, skewness, and the kurtosis (the 
latter only up to the first (T-correction) they give, 

a^ f» af + 1.27 at + 2.02 erf + O (erf) 
Ss ^ 4 + 4.69 erf + 13.93 erf + O (erf) 

Si ^ 30.22 + 98.51 erf + O (erf) . (D5) 



SSZA 


Unsmoothed 




Smoothed 






7 = 


7=-l 


7 = -2 


7 = -3 




n = —3 


n = -2 


n = -1 


n = 


S2A 


0.78 


0.47 


0.78 


1.69 


S3.Q 


4 


3 


2 


1 


S3.2 


4.37 


0.94 


-0.15 


-1.13 


Sifi 


30.22 


15.89 


6.22 


1.22 


Sa.2 


92.81 


13.85 


-0.96 


-1.82 



Table Dl. Values for the higher-order perturbative contributions 
for the SSZA for the unsmoothed (n = —3) and smoothed (n = 
—2,-1,0) density fields, for a top-hat window and a power-law 
power spectrum. 



These results are to be compared to our results from the SC 
model in the perturbative regime (truncated at the same 
order), 

o-^ ^ af + 0.78 erf + 0.76 erf + O (af) 

53 ^ 4 + 4.37 erf -I- 9.41 af + O (af) 

5*4 « 30.22 + 92.81 af -f- 331.0 af + O {af) , (D6) 

which differ in 7 % in the corrective term for Ss , 6 % in the 
corrective term for 5*4, and up to 40 % in the one- loop term 
for the variance. The two-loop contributions for 5*3 differ in 
about 30 % and just give the right order of magnitude for 
the variance. The two-loop contribution to 54 {^ 331 af) 
must be taken just as an estimate of the actual value for this 
coefficient in PT, since there are no analytic results available 
to compare to. All the above quoted differences must be 
attributed to the tidal effects as argued before what gives 
further support to the view that the hierarchical amplitudes 
are essentially of a shearless nature unlike the (unreduced) 
one-point cumulants. 

For the smoothed density field in the ZA to PT there is 
only one result available concerning the skewness for n = — 2 
(see S97), 

S3 ~ 3 + 0.82 af + 0{af), 



(D7) 



while within the SSZA dynamics, we obtain. 



5-3 ^ 3 + 0.94 af + 0(af), (D8) 

which means a 15% negative contribution from the shear 
for the first corrective term for that particular value of the 
spectral in dex. 
In Fig 



Dl 



we display the deviations from the tree-level 
values for the variance, skewness and kurtosis for the SSZA 
up to 2nd perturbative contribution (one- loop). In line with 
the arguments pointed out for the spherically symmetric 
approximation to the exact PT (the SC model), here the 
smoothing effects tend to diminish the non-linear correc- 
tions as well. We have checked that the 3rd order contribu- 
tion follows qualitatively the same behavior although non- 
linearities are typically larger. On the other hand, the lack 
of a critical index (vanishing non-linearities) for the vari- 
ance gives further support to our claim that it is due to the 
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Figure Dl. Same as Fig hi, for the SSZA. In this case, the vari- 
ance has a minimum corrective term for n Ri —1, and the hi- 
erarchical amplitudes show a vanishing first corrective term for 
n Ri —1.2. Notice that the n = 1 line is not plotted for the skew- 
ness since for that particular value the tree-level exactly vanishes. 

local nature of the SC picture and thus, to the loss of the 
previrialization effect. 
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